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Abstract 

We describe an algorithm that takes as input n points in the plane 
and a parameter e, and produces as output an embedded planar graph 
having the given points as a subset of its vertices in which the graph 
distances are a (1 + e)-approximation to the geometric distances between 
the given points. For point sets in which the Delaunay triangulation has 
sharpest angle a, our algorithm's output has O(^-n) vertices, its weight is 
0(£) times the minimum spanning tree weight where /3 = — log —. The 
algorithm's running time, if a Delaunay triangulation is given, is linear in 
the size of the output. We use this result in a similarly fast deterministic 
approximation scheme for the traveling salesperson problem. 

1 Introduction 

A spanner of a set of points in a geometric space is a sparse graph having those 
points as its vertices, and with its edge lengths equal to the geometric distance 
between the endpoints, such that the graph distance between any two points 
accurately approximates their geometric distance |12j . The dilation of a spanner 
is the smallest number 6 for which the graph distance of every pair of points 
is at most S times their geometric distance. It has long been known that very 
good spanners exist: for every constant e > and constant dimension d, it is 
possible to find a spanner for every set of n points in 0(n\ogn) time such that 
the dilation of the spanner is at most 1 + e, its weight is at most a constant 
times the weight of the minimum spanning tree, and its degree is constant pQ. 

A spanner is plane if no two of its edges (represented as planar line seg- 
ments) intersect except at their shared endpoints |17j . Plane spanners with 
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bounded dilation are known; for instance, the Delaunay triangulation is such 
a spanner [§]. However, it is not possible for these spanners to have dilation 
arbitrarily close to one. For instance, for four points at the corners of a square, 
any plane graph must avoid one of the diagonals and have dilation at least y2- 
However, the addition of Steiner points allows smaller dilation for pairs of orig- 
inal points. For instance, the plane graph formed by overlaying all possible line 
segments between pairs of input points has dilation exactly one, although its 
0(n 4 ) combinatorial complexity is high. Less trivially, in the pinwheel tiling, a 
certain aperiodic tiling of the plane, any two vertices of the tiling at geometric 
distance D from each other have graph distance D + o(D) [TB]. We define a 
plane Steiner 5 -spanner for a set of points to be a graph that contains the points 
as a subset of its vertices, is embedded with straight line edges and no crossings 
in the plane, and achieves dilation 5 for pairs of points in the original point set. 
We do not require pairs of points that are not both original to be connected by 
short paths. 

Arikati et al. [2] show how to construct a plane Steiner spanner in O(nlogn) 
time, but do not bound the total weight of the graph. Of course, spanners may 
also be constructed by forming an arrangement of line segments representing 
the edges of a nonplanar spanner graph; this planarization does not change the 
spanner's weight, but may add a large number of edges and vertices. A paper 
of Klein [TJ] on graph spanners provides an alternative basis for plane Steiner 
spanner construction. Generalizing a previous result of Althofer et al. [T] , Klein 
shows that any n-vertex planar graph with a specified subset of vertices may 
be thinned to provide a planar Steiner (1 + e)-spanner for the graph distances 
on the specified subset, with weight 0(l/e 4 ) times the weight of the minimum 
Steiner tree of the subset, in time 0((n log n) /e). Klein combined this result with 
methods from another paper [T5] to provide a polynomial time approximation 
scheme for the traveling salesperson problem in weighted planar graphs. Using 
Klein's method to reduce the weight of the geometric spanner formed by the 
arrangement of all line segments connecting pairs of a given point set would lead 
to a low weight plane (1 + e) Steiner spanner for the point set, but again with 
a large number of vertices and edges. Ideally, we would prefer plane Steiner 
spanners that not only have low weight, but also have a linear number of edges 
and vertices. 

Small and low-weight plane Steiner spanners in turn could be used with 
Klein's planar graph algorithms to derive a deterministic polynomial time ap- 
proximation scheme for the Euclidean TSP. The previous randomly shifted 
quadtree approximation scheme of Arora [3] and guillotine subdivision approx- 
imation scheme of Mitchell [H] have runtimes that are polynomial for fixed e 
but with an exponent depending on e; in contrast, Klein's method takes time 
linear in the spanner size for any fixed e. However, combining Klein's method 
with the nonlinear-size Steiner spanners described above would not improve on 
a different deterministic TSP approximation scheme announced by Rao and 
Smith 19J. Their method is based on banyans, a generalized type of spanner 
that must accurately approximate all Steiner trees, and it takes O(nlogn) time 
for any fixed e and any fixed dimension, although its details do not appear to 
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have been published yet. 

These past results raise several questions. Are banyans necessary for fast 
TSP approximation, or is it possible to make do with more vanilla forms of 
spanners? How quickly may low-weight plane Steiner spanners be constructed, 
and how quickly may the TSP be approximated? And how few vertices are 
necessary in a plane Steiner spanner? 

In this work we provide some partial answers, for planar point sets that 
are well-spaced in the sense that their Delaunay triangulation avoids angles 
sharper than a for some a. We show that, when both e and a are bounded by 
fixed constants, there exist plane Steiner (1 + e)-spanners with 0(—n) vertices 
whose weight is times the minimum spanning tree weight where (3 = 

-Mog Note that the weight depends linearly on -, improving the quartic 
dependence given by Klein's thinning procedure, which additionally has 0(n 4 ) 
vertices. Our spanners may be constructed in linear time given the Delaunay 
triangulation. In order to use our spanner for approximating Euclidean TSP, 
we may assume that our points have integer coordinates and so we can use 
the fast Delaunay triangulation algorithm of Buchin and Mulzer [TU] with fast 
integer sorting algorithms |13| ) to find the triangulation in time 0{n^/\og logn). 
Combining our spanners with the methods from Klein |15] leads to near-linear- 
time TSP approximation for the same class of point sets. 



2 Delaunay triangulations without sharp angles 

The Delaunay triangulation DT of a set S of points (called sites) is a trian- 
gulation in which the circumcircle of each triangle does not contain any sites 
in its interior. For points in general position (no four cocircular) the Delaunay 
triangulation is uniquely defined and its sharpest angle a is at least as large 
as the sharpest angle in any other triangulation. As we show in this section, 
Delaunay triangulations that do not have any triangles with sharp angles have 
two key properties: 

1. Their total weight w(DT) is small relative to the weight u>(MST) of the 
minimum spanning tree. 

2. Every point in the plane is covered by only a few circumcircles. 

We use the first property to bound the total weight of the final spanner as each 
edge we add will have length at most that of the Delaunay triangle in which it 
is embedded. In order to approximate a the length of a line, we will charge the 
error incurred to a chord given by the intersection of the line with the interior 
of a Delaunay circumcircle. Since a part of the line may be enclosed by multiple 
circumcircles, the error we charge will multiply by this factor. By bounding this 
factor, using the second property, we bound the total error. 

Lemma 1. 

1 -t- cos rv 

w(DT)< f w (a)w(MST) where f w (a) = -± . 
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Figure 1: An ordering for the edges in the proof of Lemma[T] The MST is given 
by the solid edges. The non-MST edges are dashed and, viewed in the dual of 
the plane graph defined by the DT, form a spanning tree of this dual graph. 



Proof. The proof follows closely to that of Lemma 3.1 of Klein [2]. Let T be 
the MST. (Recall T C DT.) Consider the dual graph of the plane graph defined 
by DT and refer to Figure [I] The edges DT \ T, viewed in the dual, form a 
spanning tree T* of the dual graph. Rooting T* at the vertex corresponding to 
the outside of the DT, we consider any leaf-to-root order of DT\T with respect 
to T* . Let ei, e2, . . • , efc be that ordering. 

Let H be the non-self-crossing Euler tour of T. The first non-MST edge, 
ei, makes a triangle with edges a\ and b\ of Hq. Recursively define Hi as the 
tour resulting from removing a, and 6; from iT,_i and adding e,-: 

w(Hi) = w(Hi-i) + w(ei) - w(a,i) - tu(bj) 

Since the a is the smallest angle of triangle e^a^ and e$ is longest when 
uj(oi) = w(bi), we get 

w(e.i) < (w(ai) + w(bi)) cos a. 

Combining, we get: 



w{Hi) < w{Hi- X ) + (1 - I / cos a) w{ei). 

Summing: 

w(H k ) < w(H ) + (1 - 1/ cos a) ^ wfa). 

i 

Rearranging gives: 

, n ^ w(H ) - w(H k ) 2 cos a 

> wieA < = < lo(MbT) 

^ K lJ ~ — 1 " 1 -cosa V ; 



CCS (\ 



where the last inequality follows from w(Hk) > and w(Ho) = 2w(MST). 
Since w(DT) = w(MST) + Y,i w(ei), we get the lemma. □ 
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Lemma 2. The number of Delaunay circumdisks whose interiors contain a 
given point in the plane is at most 

f e (a) = 2-K/a. (1) 

Proof. The lemma trivially holds for points that are sites. Let x be a non-site 
point in the plane. Then the Delaunay triangles whose circumcircles contain 
x are exactly the ones that get removed from the Delaunay triangulation if we 
add x to S and re-triangulate. Therefore, the number of Delaunay circumcircles 
that contain x is the same as the degree of x in the Delaunay triangulation, 
DT, of SU{x}. 

Let d be the degree of x in DT^.. Then, one of the triangles, xqr, in DT^ 
incident to x has an angle at x of at most 2n/d. Edge qr must be a side of a 
triangle qrs in DT because, after the removal of x, line segment qr is still a chord 
of the empty circle that circumscribed xqr. However since the circumcircle of 
triangle qrs contains x, this circumcircle extends at least as far from the ir-side 
of qr as the circumcircle of xqr. Therefore angle qsr is at least as sharp as angle 
qxr. So it must be that 2tt /d > a, proving the lemma. □ 




Figure 2: An illustration of the construction in the proof of Lemma [2j DT X is 
given by the dark edges and DT is given by the convex hull of DT X and the 
gray edges. The circumcircles of triangles qrx and qrs are illustrated. 
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Figure 3: Illustration of the statement and proof of Lemma [3] 

3 Portals for chords 

As we now show, it is possible to space a set of portals along an edge of a 
Delaunay triangulation in such a way that any chord of a Delaunay circumcircle 
must pass close to one of the portals, relative to the chord length. 

Lemma 3. Let AD be a chord of a circle O, let B and C be points interior to 
segment AD, and let EF be another chord of O, crossing AD between B and 
C . Then the distance from chord EF to the nearer of the two points B and C 
is at most 

\EF\ ■ \BC\ 
2min(|AB|, \CD\)' 

Proof. The points of the lemma are illustrated in Figure [3| We assume without 
loss of generality that F is on the side of AD that contains the center of O, 
as drawn in the figure; let Y be the point of O farthest from X, lying on the 
line through X and the center of O. Note that the distance from line EF to 
the closer of B and C is at most mm(|iLY|, |CX|) < \BC\/2, so it remains to 
prove that \EF\ > min(|AB|, |CD|). But if F lies on the arc between A and 
Y, then \EF\ > \FX\ > \AB\, and if F lies on the arc between Y and D then 
\EF\ > \FX\ > \CD\. In either case the result follows. □ 

Lemma 4. Let s be a line segment in the plane, and let e > 0. Then there 
exists a set P s>e of 0(- log -) points on s with the property that, for every circle 
O for which s is a chord, and for every chord t of O that crosses s, t passes 
within distance e|i| of a point in P s . e . 

Proof. Our set P S <L includes both endpoints of s and its midpoint. Refer to 
Figure |4j In the subset of s from one endpoint po to the midpoint m, we add a 
sequence of points Pi, where p\ is at distance 0(e 2 s) from p with a constant of 
proportionality to be determined later and where for each i > 1, p t is at distance 
2ed(po,Pi-i) from Pi-%. Because the distance from pg increases by a (1 + e) 
factor at each step, the set formed in this way contains 0(- log -) points. 
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If chord t crosses s between some two points pt and Pi+\ for i > 1, or between 
the last of these points and the midpoint of s, then by Lemma [3] the distance 
from t to the nearer of these two points is at most: 

1*1 • \PiPi+i\ _ 2e|t| • \poPi\ _ , , 
2\ PoPl \ " 2\ PoPl \ 61 

Otherwise, t crosses s between po andpi. Let r be the radius of O, necessarily 
at least |s|/2, and suppose that t passes within distance Sr of pq. Because of 
the choice of pi, S = 0(e 2 ). Applying the Pythagorean Theorem to the shaded 
triangle in Figure [4] we get \t\ > r\/25 — S 2 = Q,(r*/d). Combining this with 
the definition of 5 shows that t is within distance 0(VS\t\) — 0(e\t\) of p . By 
choosing the constant of proportionality for the placement of p\ appropriately 
we can ensure that this distance is at most e|t|. □ 




Figure 4: The set of points P s ^ along a chord s (horizontal line) of a circle O. 
An illustration of the final case of the proof of Lemma [4] 

We call the points in P Si£ portals. 

4 Spanning the portals within each triangle 

Within each triangle of the Delaunay triangulation, we will use a plane Steiner 
spanner that connects the portals that lie on the triangle edges. For this special 
case, we use a construction that generalizes to an arbitrary set P of points on 
the boundary of an arbitrary planar convex set K . Consider a line L and that 
makes an angle 9 with the vertical. For S £ [0, 7r/2), we say that a path is 
{9 ± S) -angle-bounded if it is piecewise linear and the smallest angle between 
each linear segment and L is at most 5. We say that a point p on the boundary 
of K is (9 ± 5) -extreme if all rays starting at p and making an angle at most S 
with L are external to K. 
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Figure 5: An illustration of the construction for Lemma [6] with 9 = 0. 



Lemma 5. Every {9 ± 5) -angle-bounded path has length 1 + 0(S 2 ) times the 
distance between its endpoints. 

Proof. The most extreme case is a path that follows two sides of an isosceles 
triangle having the endpoints of the path as base, for which the length is the 
length of the base multiplied by 1/ cos 5 = 1 + 0(5 2 ). □ 

Lemma 6. Let P be a set of n points on the boundary of a convex set K with 
perimeter I, let 9 be an angle and let 5 G [0, 7r/2). Then in time O(nlogn) 
we can construct a set S of 0(n) line segments within K , with total length 
0((£\ogn)/S), with the property that for every pointp in P there exists a (9+5)- 
angle-bounded path in S from p to a (9 + S)-extreme point of K . 

Proof. We consider the points of P in an order we will later define; for each 
such point p that is not itself (#±(5)-extreme, we extend two line segments with 
angles 9 — 5 and 9 + 5 until reaching either an extreme point of K or one of the 
previously constructed line segments. Thus, a (9 + 5)-angle-bounded path from 
p may be found by following either of these two line segments, and continuing 
to follow each line segment hit in turn by the previous line segment on the path, 
until reaching an extreme point. 

The non-extreme points of P, because K is convex, form a contiguous se- 
quence along the boundary of K. We extend segments from the two endpoints of 
this sequence, then from its median, and then finally we continue recursively in 
the two subsequences to the left and right of the median, as shown in Figure [5j 

The segments from the first two points of P contribute a total length at 
most I to S. Consider the ray s extended at an angle of, w.l.o.g., 9 + 5 from p. 
Let t be the ray extended at an angle of 9 — 5 from the most recently considered 
point p' counter-clockwise along the boundary of K from p. Consider the right 
triangle one of whose corners c is the intersection of t and s, another of whose 
corners is p and makes an angle 5 at c. (This triangle is shaded in Figure [6j) 
Let r be the right angle in this triangle. Then \pc\ = \pr\j sine) = 0(\pr\/5). 
Since a subsection of pc is added to S and the (shorter) boundary of K from p 
to p' is at least as long as pr, the length of each added segment for point p is at 
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most proportional to 1 /S times the length of the part of the boundary of K that 
extends from p to the most recently previously considered point in the same 
direction. Because of the ordering of the points, each point along the boundary 
is charged in this way for O(logn) segments, so adding this quantity over all 
points, the total length of the segments is 0((£\ogn)/S) as claimed. We may 
construct S in O(nlogn) time by using binary search to determine the endpoint 
on K of each segment. □ 




Lemma 7. Let P be a set of n points on the boundary of a convex set K with 
perimeter £, and let e > be a positive number. Then in time 0(n 2 /e) we can 
construct a plane Steiner (1 + e) -spanner for P, with all spanner edges in K , 
with 0(n 2 /e) edges and vertices, and with total length 0((£\ogn)/e). 

Proof. We choose 5 = 0(y/e) (with a constant of proportionality determined 
later). Consider the 0(1/6) angles, 9q, 9%, . . ., in [0, 2tt) that are multiples of 25. 
We apply Lemma [6] for each angle 9t with 5 as defined. We overlay the result- 
ing system of O(nfS) line segments; when two line segments from different arcs 
both have the same angle and starting point, we choose the longer of the two 
to use in the overlay. The resulting arrangement of line segments has 0(n 2 /e) 
edges and vertices and total length 0((£\ogn)/e) as required, and can be con- 
structed in time 0(n 2 /e) using standard line segment arrangement construction 
algorithms 

To see that this is a spanner, we must show that every pair (p, q) of points in 
P may be connected by a short path. Let 9 be the angle formed by the segment 
from p to q, choose i such that 9 + 5 < 9i < 6 + 36, and use Lemma [6] to find 
a (9i ± <5)-angle-bounded path pp' in the spanner from p to a (9i ± i5)-extreme 
point p' . Because of the angle bound, p' must be clockwise of q. Similarly, 
we may choose 9j within 0(5) of it + 9, and find a (9j ± (J)-angle-boundcd 
path qq' to a (9j ± £)-extreme point q' that is counterclockwise of p. These 
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Figure 7: Illustration of the proof of Lemma [7] 



two paths (depicted in Figure [7]) must cross at at least one point x, and the 
combination of the path from p to x and from x to q lies within the spanner 
and is (9 ± 0(<5))-angle-bounded. By Lemma [5j this path has length at most 
l + 0(8 2 ) times the distance between its endpoints, and by choosing the constant 
of proportionality in the definition of S appropriately we can cause this factor 
to be at most 1 + e. □ 



5 Spanner construction 

We now have all the pieces for our overall spanner construction. 

Theorem 1. Let P be a planar point set whose Delaunay triangulation is given 
and has sharpest angle a, and let e > be given. Then in time 0(n log 2 (1/ (ae))/(a 2 e 3 )) 
we can construct a plane Steiner (l+e)-spanner forP withO(n\og 2 (l/(ae))/(a 2 e 3 )) 
vertices and edges, and with total length 0(w(MST)\og(l/(ae))/(a 2 e)). 

Proof. We apply Lemma [4] to place portals along the edges of the triangulation, 
such that each chord s of a Delaunay circle passes within distance 0(ae|s|) of a 
portal on each Delaunay edge that it crosses. We then apply Lemma [7] within 
each Delaunay triangle to construct a 1 + 0(e)-spanner for the portals on the 
boundary of that triangle. 

The construction time is bounded by the time to construct the spanners 
within each triangle. Since there are 0(log(l/(ae))/(ae)) portals on each trian- 
gle, the time to construct the spanner for a single triangle is 0(log 2 (l/(ae))/ (a 2 e 3 )) 
and the total time over the whole graph is 0(nlog 2 (l/(ae))/(a 2 e 3 )). This bound 
also applies to the number of vertices and edges in the constructed spanner. By 
Lemma [TJ the total perimeter of the Delaunay triangles is 0(w(MST)/a 2 ), and 
combining this bound with the length bound of Lemma [7] gives total length 
0(w(MST) log(l/(ae))/(a 2 e)) for the spanner edges. 

To show that this is a spanner, we must find a short path between any 
two of the input points p and q. By Lemma |3j the line segment pq passes 
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within distance 0(ae|s|) of a portal on every Delaunay edge that it crosses, 
where s is the chord of one of the Delaunay circles for the crossed edge. By 
Lemma [2j the total length of all of these chords is 0(\pq\/a), so we may replace 
pq by a polygonal path that contains a portal on each crossed Delaunay edge, 
expanding the total length by a factor of at most 1 + 0((ae)/a) = 1 + 0(e). 
Then, by Lemma [7] we may replace each portal-to-portal segment in this path 
by a path within the spanner for the portals in a single Delaunay triangle, again 
expanding the total length by a factor of at most l+0(e). By choosing constants 
of proportionality appropriately, we may make the total length expansion be at 
most 1 + e. □ 



6 Approximating the TSP 

An algorithm of Klein [15] provides a linear time approximation scheme for the 
traveling salesperson problem in a planar graph. Its first step is to find a low- 
weight spanner of the graph. A subsequent paper, also by Klein [2] describes 
an algorithm that, given a planar graph G and a subset S of the nodes, finds a 
subgraph of G whose weight is 0(e~ 4 ) times that of the minimum-weight tree 
spanning S and that is a (1 + e)-spanner for the shortest-path metric on S [M] . 
This subset spanner construction can be substituted for the first step of Klein's 
approximation scheme, resulting in an algorithm for approximating the TSP 
on the subset S. However, in this more general result, the spanner construc- 
tion takes time O(nlogn), so the total time for the approximation scheme is 
O(nlogn) for any fixed e > 0. 

The first step for approximating Euclidean TSP is to round the coordinates 
of the sites to their nearest integer coordinates on a sufficiently fine grid. Doing 
so allows us to take advantage of 0(ny/\og logra) Delaunay triangulation [TU] 
made possible by fast integer sorting [15] . We may then substitute our own 
faster low-weight spanner construction for the first step of the approximation 
scheme. The remaining steps of the approximation use only the facts that the 
points we are seeking to connect into a tour are vertices in a planar graph, and 
that the whole graph has total weight proportional to the minimum spanning 
tree of the given points. Thus, we obtain the following result: 

Corollary 1. For any fixed a and e, we may find a (1 + e)- approximation to 
the optimal traveling salesman tour of sets of n points in the plane with sharpest 
Delaunay triangulation angle at most a in time 0(n) plus the time needed to 
construct the Delaunay triangulation. 

It would also be possible to design a TSP approximation scheme more di- 
rectly using a framework used by Borradaile, Klein and Mathieu [5] to solve the 
Steiner tree problem; details on how this framework applies to TSP were given 
by Borradaile, Demaine and Tazari [7] in generalizing the planar framework to 
bounded-genus graphs. Their algorithm, as interpreted for point sets in the Eu- 
clidean plane, would partition the triangles of the Delaunay triangulation into 
layers according to their depth from the infinite face in the dual graph so that the 
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sum of the lengths of line segments common to different layers is an e fraction of 
the optimal solution. This can be achieved with depth f w (a)/e = O (^); each 
layer has tree-width polynomial in this depth. The problem is then solved using 
dynamic programming, where the dynamic programs are additionally indexed 
by the portals. The base cases are made to correspond to the triangles in which 
the intersection with any tour can be enumerated. The size of the dynamic 
program is therefore bounded singly-exponentially in 1/e and 1/a. This task 
is slightly easier in the geometric setting than in the planar graph setting as 
computing shortest paths is trivial. 

7 Looking ahead 

The most obvious question posed by this work is: how do we remove the de- 
pendence on a? The dependence on a appears in two places: in the number 
of circumcircles that enclose a point and in the weight of the Delaunay trian- 
gulation. The former affects the error incurred by using portals and the latter 
affects the weight of the final spanner. We believe that it should be possible 
to remove these dependencies on a by treating groups of skinny triangles as a 
single region. In fact, using this idea, we are able to remove each dependency 
separately, but not together. By removing a long edge connecting two skinny 
triangles, we reduce the number of portals we must reroute through, but the 
number of skinny triangles that define a given region could be many, and adding 
the edges to build the spanner within this region will depend on this number. 
On the other hand, we could only consider regions defined by a small number 
of triangles, but this may not be enough to reduce the number of circumcircles 
a chord is within. 

An alternative approach to removing this dependence would be to augment 
the input to remove all sharp angles from its Delaunay triangulation, but this 
may sometimes need a number of added points that cannot be bounded by a 
function of n [5] . A construction based on quadtrees shows that every point set 
may be augmented with 0(n) points so that the Delaunay triangulation has no 
obtuse angles [5]; the resulting triangulation may also be modified to have the 
bounded circumcircle enclosure property, despite having some sharp angles, and 
may be constructed as efficiently as sorting [6] . Applying our spanner construc- 
tion method to the augmented input would allow us to completely eliminate the 
dependence on a in the time and output complexity of our spanners, but at the 
expense of losing control over their total weight. Once a spanner is constructed 
in this way, Klein's subset spanner [T3] can be used to reduce its weight, allow- 
ing it to be used in an algorithm to approximate the TSP for arbitrary planar 
point sets in time 0{n\ogn) for any fixed e > 0, but this does not improve on 
the time bound of Rao and Smith [TS] . 

Unlike in the methods of Arora [3] , Mitchell [TB] , Rao and Smith [T^] and 
Borradaile, Klein and Mathieu [5], the approximation error in our method is 
charged locally as opposed to globally. In the quad-tree based approximation 
schemes, the error incurred is charged to the dissection lines that form the quad 
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tree. In the planar approximation-scheme framework for Steiner tree, the error 
incurred is charged to an 0(MST)-weight subgraph called the mortar graph 
which acts much like the quad-tree decomposition. Our charging scheme is 
much more similar to that used by Klein for the subset tour problem in planar 
graphs [14] . However, in applying the planar approximation-scheme frameworks 
of either Klein or Borradaile, Klein and Mathieu, some error is incurred in par- 
titioning the graph into pieces of bounded treewidth. This error is proportional 
to the graph that is partitioned, which in our case is either the spanner (for 
Klein's scheme) or the triangulations (for Borradaile et. al's scheme). This er- 
ror is indirectly related to OPT by way of the O(MST) weight of the spanner 
and triangulation. By current techniques, this source of error does not seem 
avoidable. 

Finally, our spanner construction more closely ties Euclidean and planar 
distance metrics together. By unifying the approximation schemes in these two 
related metrics, it may be possible to generalize these methods to other two- 
dimensional metrics. 
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